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We report on the formation of the high-energy empty-state resonance in the electronic spectrum 
of the iron adatom on the Pt(lll) surface. By using the combination of the first-principles methods 
and the finite-temperature exact diagonalization approach, we show that the resonance is the result 
of the valence fluctuations between atomic configurations of the impurity. Our theoretical finding is 
fully confirmed by the results of the scanning tunneling microscopy measurements [M.F. Crommie 
et al., Phys. Rev. B 48, 2851 (1993)]. In contrast to the previous theoretical results obtained by 
using local spin density approximation, the paramagnetic state of the impurity in the experiment is 
naturally reproduced within our approach. This opens a new way for interpretation of STM data 
collected earlier for metallic surface nanosystems with iron impurities. 


I. INTRODUCTION 

The modern scanning tunneling microscopy (STM) 
experiment^ are aimed to reveal, identify and operate 
excitations of surface nanosystems, which is of crucial 
importance for understanding the basic phenomena in 
quantum physics such as Kondo resonance,magnetic 
anisotropy, exchange interactional between adatoms and 
for constructing novel information devices on the atomic 
level. Typically, these excitations present the smallest en¬ 
ergy scale in the system, they are of a few millivolts and 
can be reproduced by solving Heisenberg-type models .1^ 

However, the picture of the STM experiment in the 
case of an adatom on a metallic surface is not complete 
without considering high-energy resonances of hundreds 
meV.I^ They carry information concerning inter-orbital 
charge or spin excitations in the impurit^EHJ ^^d can be 
used to study electronic structure of the adatom. In this 
sense there is one important example of such high-energy 
excitations that is a pronounced peak at 0.5 eV above the 
Fermi level in the tunneling spectra of the iron adatom 
deposited on Pt(lll),IIl Pd(lll)P or W(110)P^ surfaces. 
Importantly, it was proposed that the peak can be used 
to identify iron species in surface nanosystems. 

The theoretical description of the high-energy excita¬ 
tions such as 0.5 eV peak for iron adatom on a metal¬ 
lic surface is a complex methodological and numerical 
problem. Within the Lang’s modeP^that establishes the 
connection between experimental STM spectrum and the 
local density of states at the site of the adatom, this peak 
can be attributed to the 45 states of the Fe adatom. How¬ 
ever, the ab — initio calculations based on the density 
functional theory do not confirm such a scenario. As we 
will show below the non-spin-polarized density of states 
exhibits a Sd peak in the iron adatom spectral function at 
the Fermi level (Figj^. Above Ep the density of states is 
a decreasing function without any features. The intensity 
of the 45 states is very small. 

As it was demonstrated in Refl9] the experimentally 
observed excitations can be reproduced in the framework 
of the spin-polarized local density approximation (LSDA) 


calculations. For that one should take into account the 
shift of the spin-down 3d states to higher energies above 
the Fermi level due to the ^in splitting. Then within the 
Tersoff-Hamann approaclP^ the STS is associated with 
the scattering of the s — p electrons at spin-down 
states.l^ 

Thus to describe the experiment one should define 
some ordered magnetic state of the system in the ab — 
initio calculations. As we will show below, this scenario 
is fulfilled in the case of Fe/Pt(lll). At the same time, 
the experimental conditions (low temperatures and zero 
magnetic field) correspond to the paramagnetic state of 
Fe/Pt(lll). To describe the system in this regime, a real¬ 
istic five-orbital impurity Anderson model is constructed 
and solved by using the finite-temperature exact diag- 
onalization method. The correlated density of states 
shows a peak at 0.5 eV above the Fermi level in accor¬ 
dance with the experiment .III The composition analysis of 
the eigenvectors has shown that the 0.5 eV resonance is 
originated from the valence fluctuations, which are com¬ 
bined effect of the intra-atomic exchange coupling on 
the impurity and strong hybridization with the surface 
states. We also show that the position of the resonance 
is sensitive to the coupling with other adsorbates. 


II. LDA AND LSDA RESULTS 

The first step of our investigation was to define an equi¬ 
librium atomic structure of the Fe/Pt(lll) nanosystem. 
For that we performed first-principles molecular dynam¬ 
ics simulations W using the Vienna ab initio simulation 
package (VASP)P^^^ within local density approximation 
(LDA). In these calculations the energy cutoff of 300 eV 
in the plane-wave basis construction and the energy con¬ 
vergence criteria of 10“^ eV were used. The atomic po¬ 
sitions were relaxed with residual forces less than 0.01 
eV/A. We used the PAW-PBE exchange-correlation po¬ 
tential as described in ReflTSl 

The simulations of the atomic structure of Fe adatom 
on the Pt(lll) surface were carried out within a supercell 
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FIG. 1. Schematic representation of the Fe/Pt(l 11) (left) and 
FeH/Pt(lll) (right) nanosystems simulated in this work. 


approach. The supercell contains a three-layered (4 x 4) 
Pt(lll) surface, iron atom and vacuum region of 13.4 A. 
A simplified illustration of the simulated Fe/Pt(lll) sys¬ 
tem is given by Fig{^ (left). The lattice constant for Pt 
lattice was chosen to be 3.92 A that is the experimental 
value of the lattice constant for the bulk fee PtP The 
unit cell parameters were fixed during relaxation proce¬ 
dure. The obtained vertical distance between Fe atom 
and Pt surface of 1.61 A is in good agreement with the 
reported values 

Fig gives the partial densities of states obtained by 
using the local density approximation. There is a peak at 
the Fermi level. The width of the peak is about 1 eV that 
is much larger than one would expect for a Kondo sys¬ 
tem. In contrast to the results obtained for Co/Pt(lll) 
(Refl7|) we do not observe a strong orbital polarization of 
the LDA spectra. The integration of the density of states 
gives the Fe-3d occupation of 6.5. One can see that LDA 
density of states does not reveal any feature at the en¬ 
ergy of the experimental resonance of 0.5 eV. 4s states 
of iron can be also excluded from the consideration since 
they give a negligible contribution to the spectral func¬ 
tion close to the Fermi level. As we will show below the 
obtained LDA spectra can be used to extract the hopping 
integrals between impurity and surface states. The latter 
is important to construct a realistic Anderson model for 
Fe/Pt(lll) system. 

The account of the spin-polarization within local spin 
density approximation leads to a significant change of the 
physical properties of Fe/Pt(lll). In agreement with the 


results of previous works the spin splitting produces the 
peak at 0.5 eV in the spin-down channel (FigJ^. We ob¬ 
tain the magnetic solution with moments of the 3d, As 
and Ap iron shells that are Md = 3.2 /i^, Mg = 0.04 pB 
and Mp = 0.034 respectively. The value of the total 
moment, 3.3 pb is in good agreement with the previous 
results^ and smaller than A pB one would expect for 
the isolated iron being in the atomic configuration. A 
strong hybridization of the Fe and Pt states results in a 
partial magnetization of the surface, the induced mag¬ 
netic moment of the surface can be estimated as 1.6 pb- 
We observe the peak for spin-down 3d states of iron at 
the energy of the experimental resonance. These elec¬ 
tronic and magnetic properties of Fe/Pt(lll) agree with 
that reported in the previous works 

From Figone can see that the gravity center of the 
occupied spin-up states shifts from the Fermi level to the 
energies of about -3.5 eV. They are strongly hybridized 
with the surface state and becomes completely delocal¬ 
ized in comparison with the LDA picture. At the same 
time, the width of the LSDA peak above the Fermi level 
is smaller then that in LDA, which indicates the increase 
of the localization of the spin-down 3d orbitals. The spin 
splitting that can be estimated from Figj^is in reasonable 
agreement with model estimate /^^((n^) — (ri^)), where Id 
is the Stoner parameter of about 1 eV. Interestingly, the 
number of the 3d electrons obtained in the LSDA solu¬ 
tion, = 6.13 is smaller than the LDA value of 

6.5. Since the spin-up states are almost occupied, the 
change of the total number of the electrons is related 
to the change of the hybridization of the iron spin-down 
states with the surface. 

Thus the 0.5 eV resonance observed in the STM 
experimeniP can be reproduced and explained by using 
the results of the LSDA calculations, however, for that 
one should assume the non-zero magnetization of the im¬ 
purity. In a real experimeniPl the impurity is in the para¬ 
magnetic state and below we present the results of An¬ 
derson model simulations that describe the Fe/Pt(lll) 
system without magnetic ordering. 



FIG. 2. Partial densities of states obtained from LDA calcu¬ 
lations. Blue bold and red thin lines denote 3d and As states 
of the iron adatom. The intensity of the Fe-4s is multiplied 
by ten. The dashed line corresponds to the Fermi level. The 
dotted line denotes the energy of the experimental resonance. 
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FIG. 3. Spin-resolved densities of states of the iron adatom 
obtained from LSDA calculations. The dashed line corre¬ 
sponds to the Fermi level. The dotted line denotes the energy 
of the experimental resonance. 
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III. ANDERSON MODEL 

To take into account the paramagnetic state and dy¬ 
namical electron-electron correlations in Fe/Pt(lll) we 
have constructed and solved the following Anderson 
model 

^Fe — ~ ^p^pa^pcr 

+ 2 '^ijklaa' ^ijkld^a^^cr'dlcj' dkcr 

+ ^ipcj {yipd\cr^po- + , ( 1 ) 

where and are energies of the impurity and surface 
states, d\^ {dia) and cj^(cpcr) are the creation (annihila¬ 
tion) operators for impurity and surface electrons, Vip is 
the hopping integral between impurity and surface states, 
/i is the chemical potential and Uijki is the Coulomb ma¬ 
trix element. The impurity orbital index i (j, k, 1) runs 
over the 3d-states {xy^ yz, 3z‘^ — r^, xz^ x‘^ — y‘^). A re¬ 
alistic simulation of the Fe/Pt(lll) system requires an 
accurate definition of these Hamiltonian parameters. 
Definition of the parameters. The energies and hop¬ 
pings Vip were calculated within the minimization of the 
LDA hybridization functions presented in Figby using 
the following expression, 

\Y. |2 

= ( 2 ) 
UJ Cp 
P=1 ^ 

where Np is the number of the effective orbitals describing 
the surface. The main limitation of the exact diagonal- 
ization approach is the number of the effective orbital in 
the electronic Hamiltonian. In our study we have simu¬ 
lated five 3d orbitals of the adatom. Depending on the 
symmetry each impurity orbital is connected with a cer¬ 
tain number of the surface levels. For instance, we use 
Np = 2 for xz{yz) and x‘^ — y‘^{xy) orbitals. Since the 
hybridization function of 3z‘^ — orbital demonstrates 
more complicated structure than others then for this or¬ 
bital we used Np = 3 bath states. Thus the total number 
of the effective orbitals in Eq.(l) is equal to 16, which 
corresponds to the maximal occupation of 32 electrons. 

The impurity orbital energies, were varied within the 
expression for the impurity bath Green’s function 

Gi(a;) = [w - ei - Ai(w)]“^ (3) 

to reproduce the LDA occupations for 3d states of iron 
described in the previous section. 

In turn, the elements of the Coulomb interaction ma¬ 
trix, Uijki were defined by using effective Slater integrals 
that related to the averaged on-site Coulomb interac¬ 
tion, Ud and intra-atomic exchange interaction Jh as 
described in Refl22l 

Correlated speetral funetions. The constructed Anderson 
model was solved by using finite-temperature exact di- 
agonalization solver. To analyze the excitations near the 
Fermi level we have calculated 26 lowest eigenvalues and 
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FIG. 4. The imaginary part of the hybridization functions 
(dashed line) obtained from LDA calculations in comparison 
with fitting results by using Eq. (blue solid lines). 


eigenvectors of the Anderson Hamiltonian. The main re¬ 
sult of our investigation that is the comparison of the 
correlated spectral function and experimental densities 
of states extracted from the STM spectrsPl is presented 
in Figj^ One can see that the position and width of the 
theoretical peak are in excellent agreement with the ex¬ 
perimental results. The analysis of the partial densities 
of states shows that all the 3d orbitals equally contribute 
to the peak at 0.5 eV above the Fermi level. 

We found that the variation of Ud from 2 eV to 4 
eV and Jh from 0.6 eV to 0.8 eV does not significantly 
change the correlated spectrum near the Fermi level. One 
should mention that the choice of the chemical potential, 
y is of crucial importance, since it controls the position 
of the resonance above the Fermi level. If one chooses 
the value of y in such a way to reproduce the LDA total 
number of 3d electrons that is 6.5, then the resonance 
is at 0.3 eV. The solution with correct position of the 
resonance at 0.5 eV above the Fermi level corresponds to 
6.13, which agrees with the LSD A solution. 
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FIG. 5. The spectral function of the iron impurity atom calcu¬ 
lated by using the exact diagonalization approach (blue line). 
The dotted lines correspond to the densities of states deter¬ 
mined in the STM experiment Ref{T]with different tips. 
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FIG. 6. The calculated spectrum of the Anderson model. 
Arrow denotes the excitation corresponding to the peak at 0.5 
eV above the Fermi level. No is the number of the electrons 
in the cluster, No = 27. 


To understand the microscopic origin of the peak we 
analyzed the low-energy spectrum of the Anderson model 
presented in Figj^ There are three excited states at 0.48 
eV, 0.49 eV and 0.56 eV. The excitations from the ground 
state to first excited states correspond to the experimen¬ 
tal resonance above the Fermi level. These excitations 
lead to the change of the total number of the electrons 
in the system that is increased by one. The change of 
the impurity properties at this transition can be traced 
by calculating the distribution of the atomic configura- 



FIG. 7. Weights of the impurity atomic configurations com¬ 
prising the ground and first excited states. 


tions of the impurity for each eigenstate (Fig0. For the 
ground state we obtain 0.3do + 0.47(io + 0.2do + 0.02(io, 
where the subscript index denotes the ground state. 
The iron impurity is mainly in configuration, which 
follows from the Hund’s rules for isolated atom. A 
strong coupling with the substrate leads to the valence 
fluctuation^^ to and d^ configurations within the 
ground state. In turns, the first excited state at 0.48 
eV is characterized by the following composition of the 
atomic configurations 0.56df -b O.SSdl + O.OGdf. Compar¬ 
ing the computed compositions for the ground and first 
excited states we observe the redistribution of the atomic 
configuration weights of the impurity, which corresponds 
to the do ^ di and dg ^ dl transitions. 

Experimentally, there is also a shallow rise in the den¬ 
sity of states at approximately -0.35 eV, which can not be 
reproduced by using the LSDA approach (FigJ^. This 
experimental feature can be associated with the excita¬ 
tion at -0.7 eV in the theoretical spectrum of the Ander¬ 
son model. The latter has the symmetry of the in-plane 
xy and x‘^ — y‘^ states. 

Thus we obtain two different descriptions of the exper¬ 
imental STM resonance at 0.5 eV for Fe/Pt(lll). The 
first one is the LSDA solution where the peak is the re¬ 
sult of the splitting between spin-up and spin-down 3d 
states of the adatom. On the other hand the solution of 
the five-orbital impurity Anderson model revealed a com¬ 
plex multiplet structure of the iron atom at low temper¬ 
atures, which is the result of the interplay between intra- 
atomic exchange coupling and strong impurity-surface 
hybridization. 

Importantly, one can find a connection between the 
LSDA and Anderson model results. For an isolated atom 
with a partially filled shell it was shown by Slatei!^ that 
the difference between one-electron energies with spin- 
up and spin-down is proportional to the terms that re¬ 
main even if there is no difference between spin-up and 
spin-down orbitals. Based on this result he found a con¬ 
nection between the spin-splitting energy computed in 
an one-particle approach and the difference between the 
average energy of all multiplets with S and those with 
S — 1 (for the isolated iron 5=2) calculated by using a 
many-particle method. In our case these results should 
be revised by taking into account a strong hybridization 
with the surface, which leads to the valence fluctuation 
in the Fe/Pt(lll) system. Since a further quantitative 
comparison of the LSDA and Anderson model energies 
requires the calculation of high-energy (^ 3-5 eV ) eigen¬ 
states of the Anderson Hamiltonian we left it for a future 
investigation. 

Magnetic susceptibility. To complete the picture of the 
magnetic properties of Fe/Pt(lll ), we calculated the 
spin-spin susceptibilities for 3d states of the iron adatom. 


xf(^) 


Z uj %5 Eji — Er^/ 

nn' 


( 4 ) 

where Z is the partial function, E^ is the eigenvalue of 
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the Hamiltonian, Eq.(l). One can see that the lowest 
resonance in the spectral function of x(co’) is at 1.7 eV 
(FigJ^. From FigJ^it follows that the corresponding ex¬ 
citation does not change the total number of the electrons 
in the system. On the level of the correlated density of 
states (Fi gJH it refers to the excitation below the Fermi 
level at the corresponding energy. 


IV. IRON-HYDROGEN DIMER 

In the previous work^^^^^^^ it was proposed that the 
resonance above the Fermi level can be used for detec¬ 
tion of the iron atoms in surface nanostructures. To 
check this proposition we consider the situation when 
the iron atom is hidden from the STM tip by another 
atom (Figiii right). Such a configuration can be real¬ 
ized experimentally.!^ For that the hydrogen atom was 
deposited atop the iron impurity. By using the first- 
principles molecular dynamics simulations the equilib¬ 
rium distance between H and Fe atoms was obtained to 
be 1.6 A. Our LDA calculations show that the coupling 
with hydrogen does not change the iron partial densities 
of states for the xy^ ^ xz, yz states. At the same 

time, the 3z‘^ — atomic orbital of Fe and s orbital of 
H form a molecular orbital. From Figone can see that 
the bonding and antibonding states are located at -2.4 
eV and -fO.5 eV, respectively. 

The Anderson model describing the FeH/Pt(lll) is 
given by 

'HfcH = 'Hpe + 'Hh-Fc + 'Hh, ( 5 ) 

where Hh = ^s'^sa describes s-states of hydrogen. 

sa 

'Hh-Fc corresponds to interactions between s-states of 
hydrogen and d-states of iron and can be written as fol¬ 
lows: 

'Hh-Fc = ^ (Visdl^asa + H.c.^ . ( 6 ) 

sia 

Here is energy of hydrogen state. Vis is the hopping 
between Fe and H states, aj (ag) is the creation (annihi¬ 
lation) operator for H states. 

To define the parameters of the Anderson model de¬ 
scribing the FeH/Pt(lll) system we assume that s-states 



FIG. 8. Spectral function of the total spin-spin susceptibility 
calculated by using Eq.4. 



Energy (eV) 


FIG. 9. Gomparison of the 3z‘^ — densities of states calcu¬ 
lated for Fe/Pt(lll) (blue dashed line) and FeH/Pt(lll) (red 
solid line) by using the LDA method. 



FIG. 10. Gorrelated densities of states for FeH/Pt(lll). Blue 
line corresponds to the summary density of xy, yz, xz and 
— y^ states. Red line denotes 3z‘^ — states. Dashed 
brown line denotes the spectral function of the hydrogen atom 
deposited atop the iron impurity. 


H is hybridized only with 3z^ — states of iron atom. 
Based on these assumption we minimize the hybridiza¬ 
tion functions of the 3z‘^ — orbital with two bath or¬ 
bitals for Pt and single effective orbital for hydrogen. 
The hopping integral between hydrogen and iron orbitals 
equals to Vis = 1-56 eV. 

The correlated spectral functions of the FeH/Pt(lll) 
system is presented in FigfTol One can see that there two 
resonances above the Fermi level. Similar to Fe/Pt(lll) 
the first one is at 0.5 eV which is due to the xy, x‘^ — y‘^, 
yz and xz states. The second one at 1 eV is originated 
from the iron-hydrogen hybridization. Thus, in contrast 
to Fe/Pt(lll), the resonance of the 3z‘^ — states is 
shifted to +1 eV. We observe the excitation at the same 
energy in the hydrogen spectral function that will mainly 
contribute to the STM spectrum in the simulated config¬ 
uration. 
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V. CONCLUSIONS 


In summary, we have studied the electronic structure 
and magnetic properties of the Fe/Pt(lll) system by 
means of the first-principles calculations based on the 
density functional theory and model Anderson impurity 
Hamiltonian that was solved by using exact diagonaliza- 
tion approach. It was found that both approaches suc¬ 
cessfully reproduce the main experimental feature that 
is the peak at 0.5 eV above the Fermi level in the STM 
spectrum. While within LSDA the peak originated from 
the spin splitting of the 3d shell, on the level of the An¬ 


derson model the excitation is associated with the valence 
fluctuations between atomic configurations of iron. The 
connection between these results on the basis of the semi¬ 
nal work by Slatei!^ is discussed. In addition, we predict 
the shift of the STM resonance to 1 eV in FeH/Pt(lll). 
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